Universal computation by quantum walk 
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In some of the earliest work on quantum mechanical computers, Feynman showed how to im- 
plement universal quantum computation by the dynamics of a time-independent Hamiltonian. I 
show that this remains possible even if the Hamiltonian is restricted to be a sparse matrix with all 
entries equal to or 1, i.e., the adjacency matrix of a low-degree graph. Thus quantum walk can be 
regarded as a universal computational primitive, with any desired quantum computation encoded 
entirely in some underlying graph. The main idea of the construction is to implement quantum 
gates by scattering processes. 



I. INTRODUCTION 

Quantum computers are hypothetical devices that rep- 
resent and process information according to the princi- 
ples of quantum mechanics. If built, they could solve 
certain problems, such as factoring integers dramat- 
ically faster than classical computers. 

While the earliest algorithms for quantum computers 
were based on Fourier sampling, the concept of quantum 
walk 0, Q subsequently led to a new class of quantum 
algorithms. A quantum walk is a quantum mechanical 
analog of a classical random walk. By exploiting inter- 
ference effects, quantum walks can outperform random 
walks at certain computational tasks: there are black-box 
problems for which quantum walk provides an exponen- 
tial speedup over classical computation 0, Q , and many 
quantum walk algorithms achieve polynomial speedup 
over classical computation for pro blems of practical in- 
terest @, 0, H ©, Ea HE E3, E H HE d . 

There are several ways to define a quantum analog of a 
random walk, but perhaps the simplest is the continuous- 
time quantum walk generated by the adjacency matrix 
A of some graph Q. This quantum walk takes place on 
a Hilbert space spanned by orthonormal basis states cor- 
responding to vertices of the graph, and the evolution of 
the walk for time t is described by the unitary opera- 
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In other words, the quantum walk is simply 
the Schrodinger dynamics of a particle allowed to hop 
between adjacent vertices of the graph. 

The goal of this article is to explore the power of quan- 
tum walk as a general model of computation. In partic- 
ular, I show that even a rather restricted kind of quan- 
tum walk is universal for quantum computation, meaning 
that any problem that can be solved by a general-purpose 
quantum computer can also be solved by a quantum 
walk. (This parallels the situation for adiabatic evolu- 
tion, which is also universal [13, d, EHH.) There are 
at least two motivations for demonstrating this universal- 
ity. First, it shows that quantum walk is computationally 
powerful, since in principle, any quantum algorithm can 
be recast as a quantum walk algorithm. Indeed, these 
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two models are equivalent, since the quantum walk on 
any sufficiently sparse graph can be efficiently simulated 
by a universal quantum computer d, [U [23, HU . Not 
only does this provide renewed motivation to search for 
quantum walk algorithms, but the specific construction 
may suggest new algorithmic approaches. Second, this 
construction may provide tools for quantum complexity 
theory, as discussed further below. 

A related result was shown by Feynman in one of 
the earliest papers on quantum computation, which con- 
structs a Hamiltonian that can implement an arbitrary 
quantum computation 24]. While Feynman's motiva- 
tion was to give a physically reasonable model of a quan- 
tum mechanical computing device, his result can also be 
loosely interpreted as showing the universality of quan- 
tum walk. In some sense, any time-independent Hamil- 
tonian dynamics can be viewed as a quantum walk on a 
weighted, directed graph, with weights on opposite edges 
occurring in complex conjugate pairs to satisfy Hermitic- 
ity. However, the interpretation of Hamiltonian dynam- 
ics as a quantum walk on a graph is clearest when the 
graph is sparse and the the edges are unweighted. In 
Feynman's construction, the degree of the underlying 
graph grows with the size of the computation, and the 
edges must have weights. 

This article presents an alternative Hamiltonian for 
universal quantum computation. In this construction, 
the edges are unweighted, and the graph has maximum 
degree 3 (which cannot be improved, since graphs of max- 
imum degree 2 clearly cannot be universal). The main 
idea of the construction is to represent computational ba- 
sis states by quantum wires, and to implement quantum 
gates by scattering off widgets attached to (and connect- 
ing) the wires. For a circuit on n quantum bits, the graph 
consists of 2" wires, together with the attached widgets. 
Each wire can be idealized as infinitely long, but in prac- 
tice can be well approximated using only poly(n) vertices. 
The computation begins at a single vertex on the far left 
side of one of the quantum wires, and the outcome of the 
computation corresponds to the wire on the far right side 
to which the state evolves under the quantum walk. 

The remainder of this article is organized as follows. 
In Sec. HU we briefly review scattering theory on graphs, 
the main technical tool used in the paper. In Sec. IIII1 
we describe the set of universal quantum gates used in 
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the construction, and explain how to implement them 
by scattering processes at a certain fixed momentum. In 
Sec. IIV| we show how undesirable momentum compo- 
nents can be removed by a filtering process, also achieved 
by scattering. In Sec.fVl we describe how the various ele- 
ments of the construction behave under composition, and 
in Sec. IVI1 we argue that bound states can be neglected. 
In Sec. IVII1 we summarize the construction of a univer- 
sal quantum computer. Finally, in Sec. I Villi we conclude 
with a brief discussion of the results and some directions 
for future work. 



II. SCATTERING ON GRAPHS 

First, consider an infinite line of vertices, each corre- 
sponding to a computational basis state \x) with ieZ, 
where vertex x is connected to vertices x ± 1 . The eigen- 
states of the adjacency matrix of this graph, parameter- 
ized by k £ [— 7r,7r), are the momentum states \k) with 



(x\k) 



(1) 



for each igZ (normalized so that (k\k'} = 2tt 5(k — k')), 
with eigenvalues 2 cos k. 

Now consider an arbitrary finite graph G, and create 
an infinite graph with adjacency matrix H by attaching 
a semi-infinite line to each of N of its vertices. Label the 
basis states for vertices on the jth semi-infinite line as 
\x,j), with x = at the vertex in the original graph and 
x = 1,2,... moving out along the line. On each semi- 
infinite line, an eigenstate of the adjacency matrix must 
be a linear combination of states of the form ([1]) with 
momenta ±k, corresponding to an eigenvalue 2 cos k; or 
possibly of the same form but with k = in or k = in + n 
for some k > 0, corresponding to an eigenvalue 2 cosh n 
or —2 cosh k, respectively. For each j £ {1, . . . , N} and 
each k £ [— ir, 0], 1 there is an incoming scattering state 
of momentum k, denoted |/c, scT^) , of the form 



(x,j\k,scj*) 



-ikx 



R 3 (k)e 



ikx 



(x,/|fc,sc^)=T J , / (fc)e ite , jVj 



(2) 
(3) 



on the semi-infinite lines. The reflection coefficient Rj (fc), 
the transmission coefficients Tj,j/(k), and the form of 
|fc,scj*) on the vertices of G are determined by the con- 
dition H\k,scJ*) = 2cos/c|fc,sap). For any fixed fc, these 
coefficients can be found by solving \G\ linear equations. 
Together with bound states \k, bd ) (for certain discrete 
values 2 of k > that can be obtained by solving \G\ — 1 



1 Since the off-diagonal elements of the adjacency matrix are pos- 
itive, whereas in a discrete approximation to the kinetic term 
— d 2 /da; 2 they are negative, we use the convention that incom- 
ing states have negative momentum. 

2 Since \\H || is at most the maximum degree d of the graph, k < 
cosh" 1 ^). 



linear equations and one transcendental equation) of the 
form 



(x,j|«,bd ± ) =Bf(n) (±e- K ) a 



(4) 



the states |fc,sc~*) form a complete, orthogonal set of 
eigenfunctions of H that are useful for the analysis of 
scattering off the graph G. 

To calculate the propagator for scattering through G, 
we can expand the evolution in the basis of incoming 
scattering states and bound states. For evolution from 
vertex x on line j to vertex y on line j' ^ j, we have 



(yj'\e- im \x,j) 



E 

■7=1 



-2it cos k 



(y, f\k, sO (fc, sc^ \x, j) dk 



^ ' gT 2 if cos 

K,± 




hK ( 2 /,/|K,bd ± )(«,bd ± |x,j) (5) 



-2it cos k 



R*e 



-ikx \ 



+(e 



-iky 



R r e lkv )T*, t 



j' ,3 



dk 



V e T2it cosh K B± ( K )Bf (K)*(±e- K ) x+y 



K,± 



(G) 



-2it cos k 



%T^ r e lk[x+y) + T* je - ik ^ + ^)dk 



e T2it cosh K£± ( ±e -K- 



x+y 



(7) 



K,± 



where dk := dk/2ir, and where in the last step we have 
used unitarity of the S'-matrix. 3 (We often suppress the 
dependence of the transmission and reflection coefficients 
on momentum for notational convenience.) 

We will argue later that the contribution from the 
bound states can be neglected. According to the method 
of stationary phase (see for example 0, Sec. II.3]), the 
integral over k is dominated by those values where the 
derivative of the phase of the integrand vanishes. The 
second term in the integrand can be shown not to have 



Similarly to the incoming scattering states, we can define outgo- 
ing scattering states sc*~) of the form 



(x, j\k, sc*"> 



+ Rj(k)* e~ ik 



(ir.j'ifc.scj- > =T^,(fc)* e - ite , jVJ 

on the semi-infinite lines. The incoming states are related to 
the outgoing states by a transformation known as the scattering 
matrix, or S-matrix, with \k, sc~*) = Sjj/(k)\k,sc*7j). It 
can be shown that the S'-matrix is unitary, and has the form 
Sjj{k) = Rj(k) and Sj j/(k) = Tj ji(k) for j ^ j'. In particular, 
the orthogonality of columns j and j' of S implies that Tj j/Rj + 
Rj/ R*y . + 12j<t{j,j>} T M ,T j,j = °> § ivm g the cancellation in <(Tj - 
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any stationary points, and the phase of the first term is 
given by k(x + y) + argTjji(k) — 2icosfc, which is sta- 
tionary for 



x + y + £j,f(k) = v(k)t, 



where 



v(k) :— — 2cosfc 
dk 



-2 sin k 



is the group velocity at momentum k, and 



(8) 



(9) 



(10) 



is the effective length of the path through G from line jf 
to line j'. 4 Then for large x + y we have [25], Eq. 3.2] 



^/27T|c(fe*)| ' 



where k = k* satisfies (JSJ), and 



c(fc) := 2icosfc + — - argTjj'(fc). 



(11) 



(12) 



While semi-infinite lines are convenient for the pur- 
pose of analysis, they can be replaced by long but finite 
lines to give a construction based on a finite graph (cf. 
H). This replacement does not significantly change the 
dynamics since the quantum walk on a line has a maxi- 
mum propagation speed. To see this, note that in ([9]). a 
maximum group velocity of 2 is obtained at k = —tt/2. 
Alternatively, consider the propagator on an infinite line 
with adjacency matrix H: 



(y\e- im \x) = 



Jk(y-x)-2itcosk^j i , Qg\ 



= (-i)«-J„_ ie (2t) ) (14) 

where J v [t) is a Bessel function of order v. Since J v [t) 
decays exponentially in v when v = t(l + e) for any fixed 
e > 0, HH) describes a wavefront moving with speed 2. 
Thus, provided the lengths of all the attached lines are 
large compared to twice the total evolution time, the ef- 
fect of truncating the lines is negligible. 



III. UNIVERSAL GATE SET 

We now show how to implement a universal set of 
quantum gates by scattering on graphs. We use a univer- 
sal gate set consisting of the controlled-not gate together 



4 If the graph G is simply a line of I edges, then the transmis- 
sion coefficient is T(k) = e lke , and the effective length is pre- 
cisely £. In general, however, the effective length is momentum- 
dependent, i.e., the propagation is dispersive. 
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FIG. 1: Widgets used to construct a universal quantum com- 
puter. Open circles indicate vertices where previous or suc- 
cessive widgets can be attached, (a) Controlled-note gate, (b) 
Phase shift, (c) Basis-changing gate, (d) Momentum filter, 
(e) Momentum separator. 



with two single-qubit gates that generate a dense subset 
ofSU(2). 

The controlled-not gate is trivial to implement. This 
two-qubit gate exchanges the computational basis states 
1 10) and 1 11), while leaving the other two states un- 
changed. This transformation can be effected by sim- 
ply exchanging the appropriate wires, using the widget 
shown in Fig.[T]Ja). An incoming wave of any momentum 
k is transmitted perfectly through this widget, accumu- 
lating a phase of e lk . 

To implement a phase gate, we would like to apply 
some nontrivial phase to the |1) wire, while leaving the 
|0) wire unchanged. This can be accomplished by insert- 
ing the widget shown in Fig. [IJb) into the |1) wire. To 
understand this widget, consider attaching semi- infinite 
lines to its terminals, and calculate the transmission co- 
efficient for a wave of momentum k incident on the input 
terminal. We find 



T. 



(b) 



icos 2k esc 3 k sec k ' 



(15) 



whose magnitude squared is plotted in Fig. [2] In partic- 
ular, this widget has perfect transmission at k = — 7r/4, 
where T^{-i:/A) = 1 and £( b )(-7r/4) = 1. Relative to 
the effect of a straight wire of length 1, the widget ef- 
fectively introduces a phase of e 17r / 4 at this momentum. 
Combining the widget on the |1) wire with a straight wire 
for the |0) state, we see that for momenta near — 7r/4, the 
widget implements the phase gate 



U h := 



1 

e i7r / 4 



(16) 



Note that momenta far from — ir/4 (and — 37r/4) will not 
only be transmitted with a different phase, but will also 
include a substantial reflected component. However, we 
will see that the computation can be performed entirely 
with wave packets consisting of momenta near — 7r/4. 

To implement a basis-changing single-qubit gate, we 
must design a widget that includes interactions between 
different quantum wires. Such a widget is shown in 
Fig. [ljc). To characterize this widget, we calculate the 
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k 

FIG. 2: Transmission probability for the phase shift widget 
(Fig- Kb)). 




FIG. 3: Transmission probabilities for the basis-changing gate 
widget (Fig. Q^c)) with input at |0i n ) and outputs at |0 ou t) 
(solid line), |l out ) (dashed line), and |lj n ) (dot-dashed line). 



reflection and transmission coefficients for a wave of mo- 
mentum k incident on one terminal (say, the one labeled 
|0in) in Fig. Hfc); the others are related by symmetry). 
We find 



r (c) 



T, 



(c) 



(c) 

in 



T, 



(c) 



Oi„,l in 



e lfe (cosfc + isin3/c) 
2 cos k + i(sin 3fc — sin k) 
1 

2 cos k + i(sin 3fc — sin k) 

e lk cos 2k 
2 cos k + i(sin 3fc — sin k) 



(17) 
(18) 
(19) 



The corresponding transmission probabilities are shown 
in Fig. [3] At k = — n/4, the input amplitude is trans- 
formed into an equal superposition of output amplitudes, 
with no amplitude reflected back to the input chan- 
nels. The effective lengths for forward transmission are 

^o-io out (- 7r / 4 ) = ^iioutC-""/ 4 ) = 2 > so that the wid- 
get effectively lengthens the wires involved by two units. 
Considering the phases of the transmission coefficients, 
we see that transmission through the widget effectively 
performs the unitary transformation 



i 1 
1 i 



(20) 



It is straightforward to show that this gate, together with 
the phase gate (fTo) , generate a dense subset of SU(2) — 
for example, because U^U C U£ is the Hadamard gate (up 
to a global phase) [26[ ■ 

So far, we have only described how these gates act on 
one or two qubits at a time, but it is straightforward to 
embed them in a graph representing a computation on n 
qubits. For the controlled-not gate, we simply include its 
widget 2 n ~ 2 times, once for every possible setting of the 
n — 2 qubits not involved in the gate. Similarly, for the 
single-qubit gates, we include their widgets 2™ _1 times. 
As an example, Fig.Ushows the graph corresponding to a 
simple two-qubit quantum circuit. Notice that, although 
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FIG. 4: Graph implementing a Hadamard gate on the second 
qubit followed by a controlled-not gate with the second qubit 
as the control. 



the graph corresponding to an n-qubit circuit is exponen- 
tially large in n (as it must be to represent an exponential 
number of basis states), it has a succinct description in 
terms of the original circuit being simulated. 

Using only the three gate widgets (a), (b), and (c), 
we can already construct a universal quantum computer, 
provided the input state is chosen appropriately. Since 
there is no reflection at k = — 7r/4, the transmission co- 
efficients at this momentum compose multiplicatively, so 
the concatenation of gate widgets can describe an arbi- 
trary quantum circuit. If the input state is prepared in a 
narrow wave packet consisting only of momenta close to 
k = — 7r/4, the propagation of this wave packet through 
the widgets implements that circuit. However, we will 
see next that it is possible to use a much simpler start- 
ing state, corresponding to one particular vertex of the 
graph. 



IV. MOMENTUM FILTERING 

To construct a Hamiltonian that works with a simple 
starting state, we design a filter that only allows mo- 
menta near k = — 7r/4 to pass. The basic building block 
of this filter is shown in Fig. [IJd) . Unlike the widgets for 
implementing gates, this widget includes a semi-infinite 




FIG. 5: Transmission probabilities for the filter widget 
(Fig. [Hd)) with input at |in), and outputs at |out) (solid 
line) and the semi-infinite line exiting upward (dashed line). 



line, which allows undesired momentum components to 
be carried away. 

The transmission probabilities for this widget are 
shown in Fig. [51 Note that at k = — 7r/4 and k — —3ir/4, 
all amplitude is transmitted forward, whereas at other 
momenta, some of the amplitude is transmitted upward 
and some is reflected. 

To filter out all but an arbitrarily narrow range of mo- 
menta, we repeat this widget many times in series. The 
scattering properties of such a graph can be analyzed us- 
ing a transfer matrix technique described in [2] . Consider 
an infinite line of vertices with additional edges attached 
to vertices 1, . . . , ma, and an arbitrary graph (here, the 
filter widget) attached above each of those vertices. Then 
the amplitudes of the scattering state \k, sc^) satisfy 



FIG. 6: Absolute values of eigenvalues of the transfer matrix. 



[x + l\k,sc£)\ M ( (x|fc,sc£) 
(dfc,sO / Was- life, scr 



where 5 



with 



M 



y(k) 



2 cos k — y(k) — 1 
10 



(vertex above wire|fc, sc[~^) 
(vertex on wirejfc, sCjT*) 



(21) 



(22) 



(23) 



The transformation from the input amplitude to the out- 
put amplitude of the chain of md filter widgets is de- 
scribed by the matrix M md . One can show that for 



M r ' 



a b 
c d 



(24) 



5 This matrix differs from that in 0] since we are using the adja- 
cency matrix rather than the Laplacian as the Hamiltonian. 



the transmission coefficient is 

2ie -ifcm d sinfc 



(d) 



-ik 



- b + c + de ik ' 



For the filter widget, a calculation shows 6 



y(k) = 



2ik 



cos 2k 
sin k 



(25) 



(26) 



The absolute values of the eigenvalues of the correspond- 
ing transfer matrix (|22|) are shown in Fig.[6j Except when 
k is close to — tt/A or —3ir/4, one of these eigenvalues 
is bounded above 1; this eigenvalue gives the dominant 
contribution to (|25|). and the transmission coefficient is 
exponentially small in m^. However, y{— 7r/4) = 0, 
so there is perfect transmission at k = — tt/4. Since 

4nout( — 7r /^) = 2, each instance of the widget effectively 
lengthens the wire by two units at this momentum. 

Unfortunately, this filter transmits undesired momenta 
near k = — 37r/4 as well as the desirable momenta near 
k = — tt/4. Generically, two distinct momentum com- 
ponents propagate at different speeds, and hence can be 
isolated temporally. However, these particular momen- 
tum components have the same group velocity ([5]). 

To isolate the two components, we can use the widget 
shown in Fig.QJe), which has a different effective length 
for k — —it/A and k = —3ir/4. 7 The transmission coeffi- 
cient for this widget is 



(e) _ 



i(cos/c -I- cos3fc) 



sin k + 2 sin 2k + sin 3fc — sin 5fc 



(27) 



6 While y(k) is real for any finite widget, the presence of a semi- 
infinite line allows it to be complex. 

7 The symmetry of the transmission and reflection coefficients un- 
der k — > —it — k for widgets (a)-(d) holds because those graphs 
are all bipartite. Notice that widget (e) breaks this symmetry. 
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FIG. 7: Transmission probability for the momentum separa- 
tor widget (Fig. [l^e)). 



ward, respectively, as 



i? 7 



jinj; n 



3=3 
3 * 3' 

3=3' 
3 * f 



(28) 
(29) 



f° r Jj j' S {0, . . . , 2™ — l}. 8 If we place two widgets in 
series, we can compute the transmission and reflection 
matrices for the composed widget as follows: 

T 12 =T l {l-n 2 'Riy 1 T 2 (30) 

fti2 =n 1 + T 1 (i-n 2 ii 1 )- 1 Ti 2 f 1 (3i) 
r 12 =r 2 (i-^ 2 )- 1 r 1 (32) 
n 12 = n 2 + T 2 (i-n 1 n 2 )- 1 '}z 1 T 2 . (33) 



as pictured in Fig. [7J There is perfect transmission at 
k = — 7r/4 and — 37r/4 (as well as k — —tt/2). Further- 
more, the derivative of the phase of the transmission co- 
efficient gives 4n ) ou t(- 7 <7 4 ) = 4 ( 3 ~ 2 \/2) ~ °- 686 at 
k = -tt/4, and ^nloutM 71 "/ 4 ) = 4 ( 3 + 2 ^2) « 23.3 at 
fe = — 37r/4. Since the effective length of the widget is 
different for these two momenta, it serves to temporally 
separate them. 



V. COMPOSING WIDGETS 



It is straightforward to combine these widgets to sim- 
ulate an arbitrary m-gate quantum circuit: simply con- 
nect the widgets in sequence corresponding to the gates 
in the original circuit. For momentum k = —tt/4, the 
transmission coefficients of the resulting graph exactly 
implement the desired circuit. Furthermore, the effective 
lengths (fTU)) and curvatures (|12[) at k = —tt/4 simply 
add, so both are proportional to m for the overall graph. 
However, a propagating wave packet comprises a range of 
momenta, so we need to determine how close k must be 
to — 7r/4 such that the transmission gives a good approx- 
imation to the desired circuit. In turn, this determines 
how many filter widgets to include. 

To understand how gate widgets behave under com- 
position, it is helpful to view the scattering problem in 
terms of 2™ input channels scattering into 2™ output 
channels. Define matrices T, 1Z, T, TZ, describing forward 
transmission, reflection from forward to backward, back- 
ward transmission, and reflection from backward to for- 



These expressions can be obtained by constructing 
the scattering states for the composed widget out 
of scattering states for the individual widgets. For 
example, to construct |fc,sc~^)i2, we begin by con- 
catenating |fc,sc~^)i (with its output wires removed) 
and ^2iji(T\)j. m ji t \k,scp) 2 (with its input wires re- 
moved). This state fails to satisfy the eigenvalue con- 
dition only at the boundary between the widgets, so 
we attempt to correct this by including a reflection off 
widget 2: we add J2j'( T i)j^,j' ov , t (( R ^)j' oat ^jL) 1 + 
J2j"^j'( T 2)j' out ,j' c ; ut \k,sc^}i) to the state of widget 1. 
Continuing to add terms corresponding to all the possi- 
ble ways for waves to reflect between the two widgets, 
we obtain an infinite series that satisfies the eigenvalue 
condition even at the widget boundary. Summing this 
geometric series, we obtain (|517)) and (|3"Tj) , and a similar 
construction of the states |fc,sc~^ t )i2 gives ([3"2"| and l[3"3"]). 

Now we can see that the composition of two widgets, 
each with little reflection, also has little reflection. In 
particular, suppose \\7Zi\\ < #i and ||72-2||, H^H < 

S 2 . Then we have 

\\n 12 \\ < ||7ei|| + ||(i - n^y'W \\n 2 \\ (34) 

< Sx + (1 + 5182)62, (35) 

as well as a similar bound for ||72-i2||- Furthermore, 
forward transmission through the compound widget is 
nearly described by the product of the two forward trans- 



For widgets symmetric under interchanging the roles of input and 
output, as is the case for all our widgets, TZ. = TZ and T = T, 
However, this does not hold when composing distinct widgets. 
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mission matrices, since 



IT12-T1T2II = 



71 



1 



< 



1 - 7^1 

TZ 2 iZi 



n 2 ni 
- 1 



1 r 2 



1 - ^2^1 

< <JkS 2 (1 + <JkS 2 ) 



(36) 

(37) 

(38) 
(39) 

We can apply these bounds to study the transmission 
through to gate widgets for momenta near k = — tt/4. In 
particular, suppose \k + tt/4\ = 0(l/m 2 ); then ||72.|| = 
0(l/m 2 ) for each of the gate widgets. Applying the 
bounds recursively, we find that the collection of all m 
gate widgets has ||72.|| = 0(l/m) for such momenta; in 
other words, the transmission is nearly perfect. 

To filter out all undesired momenta, suppose we pre- 
cede the gate widgets by = log 9 (to 2 ) filter wid- 
gets. Then the output of the filter has exponentially 
small amplitude except for momentum components k 
with \k + tt/4| = 0(1/to 2 ) and \k + 3tt/4| = (9(1/to 2 ). 



VI. BOUND STATES 

It remains to show that bound states of the graph can 
be neglected. Since bound states decay as e~ KX , we can 
effectively ignore bound states with, say, k = fl(l/m 4 ) 
by starting the walk a distance x = 0(to 4 ) from the first 
widget. This choice increases the required running time 
of the simulation, but only by a polynomial factor. 

Strictly speaking, the above argument leaves open the 
possibility that a large number of very weakly bound 
states could influence the scattering behavior. However, 
since the bound states in question have k = 0(l/m 4 ), 
they have the very similar energies ±2 cosh k — ±(2 + 
0(l/m 8 )). Thus, for t = 0(m 4 ), the phases of the weakly 
bound states of type + in are all the same up to 
0(l/m 4 ), and similarly for the weakly bound states of 
type -. 

To have nearly identical phases for both types of bound 
states, we can simply choose an evolution time t for 



+2it 



-2ii| 



is zero. For 



which the phase difference |e 
propagation precisely at momentum k = — tt/4, corre- 
sponding to the group velocity v(— tt/4) = y2, we would 
ideally choose t — (x + £)/y/2, where i is the total ef- 
fective length of all widgets at k = —tt/4. However, 
we can vary t slightly without significantly affecting the 
contribution from the scattering states. Specifically, we 
choose t = TT[(x + e)/V2Tr\ = ((x + e)/V2)(l + 0(l/m 4 )). 
This ensures that e T2itcoshK = + 0(1/to 4 ), where the 
phase <fi is independent of k and the choice of ±. In 
other words, all weakly bound states enter with approx- 
imately the same phase. Moreover, the momentum for 
which the phase is stationary (i.e., satisfying l[5])) remains 
k* = — 7t/4 + 0(1/to 4 ), so the contribution from the scat- 
tering states is essentially unchanged. 



Now observe that at t = 0, the amplitude ([7]) is 0, and 
since the initial contribution from the scattering states 
is negligible, the initial contribution from the bound 
states is also negligible. Since the relative phases be- 
tween weakly bound states at a time t — 0{m ) can only 
change by 0(l/m ), the amplitude deviation caused by 
neglecting the bound states in {7J is 0(1 /to 4 ), which will 
turn out to be dominated by the contribution from the 
scattering states. 



VII. UNIVERSAL COMPUTER 

Overall, an arbitrary m-gate quantum circuit describ- 
ing a unitary transformation U onn qubits can be imple- 
mented by quantum walk as follows. On the input wire 
Oin, we place toj = logO(TO 2 ) filter widgets, followed by 
a momentum separation widget. We then add widgets 
for the to gates in the circuit, each with 2™ inputs and 
2™ outputs. The 2™ input wires, 2" output wires, and TOd 
filter wires are all truncated (as discussed at the end of 
Sec.|H| at length greater than 2t ~ V2(x + £) (say, length 
2(x + £)), where I is the total effective length of the fil- 
ter, momentum separation, and gate widgets. Initially, 
the computer is prepared in the state \x, 0j n ) correspond- 
ing to vertex x on input line 0, where x — 0(m 4 ) as 
discussed above. We evolve this state with the Hamilto- 
nian given by the adjacency matrix of the graph for time 
t = TT [{x + 1)/V2tt\ = 0(to 4 ). Finally, we measure m 
the vertex basis. If the outcome is on some output wire 
s € {0,1}™ (which happens with probability fl(l/m 4 )), 
we output that s; otherwise we discard the result and 
start over. 

Conditioned on obtaining a valid output, the statistics 
of this simulation closely reproduce those of the origi- 
nal quantum circuit. According to (J7]), the amplitude to 
propagate to vertex on output line s is well approxi- 
mated by 



(0,s out |e- lWt |x,O in ) » / e 



Akx — 2ii cos k 



which by (JTTJ) (and the construction of the graph) is ap- 
proximately (s|J7|0) x £1(1/to 2 ). 



VIII. DISCUSSION 

The construction described in this article shows that 
quantum walk is a universal computational primitive, in 
the sense that any quantum computation can be effi- 
ciently simulated by a quantum walk on a sparse, un- 
weighted graph. In addition to establishing the univer- 
sality of quantum walk, it is possible that ideas from this 
construction could be applied elsewhere in quantum in- 
formation processing. We conclude by briefly mentioning 
some of these possibilities. 
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One potential application is the design of quantum al- 
gorithms based on scattering on graphs. In such algo- 
rithms, we need not think of scattering as directly im- 
plementing quantum gates on one or two qubits, but 
rather, as a generic way of implementing transformations 
on quantum states. An early proposal along these lines 
was made in Q in the context of decision trees, and an al- 
gorithm for evaluating balanced binary game trees using 
similar ideas was given in (T3 | . Indeed, generalizations of 
this approach have led to quantum algorithms for evalu- 
ating broad classes of formulas [HI, [l6[ • 

Another possible application is to quantum complexity 
theory. Feynman's original quantum computer Hamil- 
tonian has been used to construct QMA-complete 
[13| and BQP-complete [H, HI problems, in addi- 
tion to showing the universality of adiabatic evolution 
Qjij]. The quantum walk construction could po- 
tentially be applied to construct new complete problems 
or adiabatic quantum computers with desirable proper- 
ties. 

Finally, this construction might be useful in the de- 
velopment of new architectures for quantum computers. 
Although we have not concerned ourselves with the phys- 
ical representation of the computer, the graph inherits a 



kind of tensor product structure from the underlying cir- 
cuit, so it might be possible to encode the system such 
that the Hamiltonian involves only local interactions be- 
tween qubits. Such a construction could potentially be 
used to design a quan tum computer that does not require 
dynamic control [H, M, M, \M, M, H H ^ . 
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